Lattice-Induced Double- Valley Degeneracy Lifting in Graphene by a Magnetic Field 
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We show that the recently discovered double-valley splitting of the low-lying Landau level(s) 
in the Quantum Hall Effect in graphene can be explained as perturbative orbital interaction of 
intra- and inter-valley microscopic orbital currents with a magnetic field. This effect is provided 
by the translational-non-invariant terms corresponding to graphene's crystallographic honeycomb 
symmetry but do not exist in the relativistic theory of massless Dirac Fermions in Quantum Elec- 
trodynamics. We discuss recent data in view of these results. 
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The recent discovery of massless charge carriers with 
linear conic spectrum - Dirac fermions (DF) in both 
graphite fl| and graphene 0, Q has prompted researchers 
to revisit many basic ideas in Solid State Physics based 
on relativistic particle physics on a lattice. The existence 
of Dirac fermions has been confirmed by direct ARPES 
0, 0| and STS measurements and by Quantum Hall 
Effect (QHE) measurements in both graphite 0, S| and 
graphene 0, S] , It has been recognized that the analogy 
between DF in graphene and relativistic massless DF in 
Quantum Electrodynamics (QED) can be used fruitfully 
to explore the properties of graphene-based systems 9] . 
One of the most important consequences of this analogy 
is the peculiar quantization of relativistic Landau Levels 
(LLs) in a magnetic field that are symmetric for positive 
(electron) and negative (hole) energies. The LL energies 
are proportional ±^/n as a function of the level number n 
away from the zero LL {n — 0) that is positioned exactly 
at the Dirac point where electron and hole spectra touch 
at zero magnetic field, H — 0. In addition to the infinite 
Landau degeneracy, each level is double spin-degenerate, 
which is a direct consequence of the relativistic (Lorenz) 
invariance of the QED equations. 

The QED double-spin degeneracy corresponds to the 
double-valley degeneracy of LLs in graphene and, to- 
gether with conventional spin-degeneracy (not considered 
by the QED analogy), produces the conductance steps of 
4e^//i [l3 observed in the semi-integer QHE 0, dou- 
ble the size of standard steps. Importantly, therefore, the 
recently discovered double-valley splitting for (at least) 
the zero Landau level [111, [13] indicates a breakdown of 
relativistic invariance in graphene. Several mechanisms 
[13] based on spontaneous symmetry breakin g d rive n by 
either long-range Coulomb interaction 14, 15l.ll6l . 17, isj . 
field-enhanced electron-phonon interaction [l9j. disorder 
13, 21, 13] or edge effects 12, 23, 24] have been proposed 



to explain this phenomenon. 

In this paper, we demonstrate that the valley gap 
opening for low-lying LLs is the intrinsic property of 
graphene-like systems. These systems have a honey- 



comb crystallographic group that is different from the 
relativistic Lorenz group in QED albeit resulting in a 
similar Dirac-like equation for non-interacting fermions 
in zero magnetic field, H = 0. The difference becomes 
apparent in an applied magnetic field when the additional 
translational-non-invariant terms accounting for interac- 
tion of microscopic intra- and inter- valley orbital currents 
with the magnetic field appear in the graphene Hamilto- 
nian. The effect of the double-valley LL splitting has, 
therefore, a much more natural explanation as a pertur- 
bative non-critical orbital splitting that is of the same 
order as the standard Zeeman spin-splitting. 

We first consider the origin and symmetry properties 
of the Hamiltonian, the spectrum, and the wave func- 
tions of conducting electrons (holes) in the vicinity of 
two crystallographically nonequivalent opposite corners 
Ki,2 (also denoted as K and K') of the hexagonal Bril- 
louin Zone of graphene at zero field, H = 0. The wave 
functions of the zero-energy states are located exactly at 

1.2 and can be linearly expanded over a 4-component 
Bloch basis (irreducible representation) of the K-point 



(1) 



(for symmetry reasons our set {^^i} is different from the 
commonly used {^^^ , ^-f ^ , ^-^^ , ^-fj) . 

It is the transformation properties of spinor-like func- 
tion vp under the action of the graphene crystallographic 



TABLE I: Transformation properties of the Bloch spinor 'I' 
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(Table |T| that define all the physical properties of charge 
carriers in graphene. Here, T1.2 are the lattice transla- 
tions; other notation are the same as in |25|. The phys- 
ical properties can be obtained either directly from the 
standard Tables of Irreducible Representations of Crys- 
tallographic Groups [26*1 or from the explicit form of ^ in 
a tight-binding approximation of the carbon pz orbitals 
marked as 7r(r) (see also Fig. [ij: 
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FIG. 1: Tfi-point Bloch functions vf^-^ and The K2- 

point Bloch functions ^f^^ and ^^'^ obtained from them 
by applying complex conjugation. 



where sa(b) = ■^-ftTi 2 — =^1 and t^^,^^ are the A 

(B) sublattice coordinates. 

Wave functions of states deviating from Ki^2 by a 
small vector k — {kx,ky) can also be expanded over 
the basis '^i{r), but with slowly space- varying envelopes 
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Ki,K2 



The energy spectrum £'(k) and the corresponding enve- 
lope functions F^{r) are the eigenvalues and eigenfunc- 
tions of the usual If fc-perturbation equation: 



H F^{r) = E{k) F^{r) 
where the iffc-perturbation Hamiltonian, 



(5) 
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(with k — — i?iV) was obtained as a most general 4x4 
matrix that is linear in k and conserves the form (^'i/^') 
under the action of the group G. 

The Hamiltonian ^ has the structure of the relativis- 
tic Dirac Hamiltonian for massless fermions with a linear 
conical spectrum (Fig. [2^): 



E{k) = ±v\k\, 



(7) 



and the corresponding system of eigenfunctions F^{r) 
that is a linear superposition (with arbitrary complex 
constants ci, C2) of two- valley plane- wave functions: 

F'^ir) = ci{±l, e'", 0, 0}e*"- + C2{0, 0, ±1, e^^}e*^ (8) 

where = avctan [k^ / ky) and the ± sign corresponds to 
the upper (lower) branch of the conical spectrum ([7] ) . 



FIG. 2; (a) Two-valley degenerate Dirac-like spectrum E{k) 
of charge carriers in zero magnetic field, H = 0. (b) Landau 
Level (LL) quantization of Dirac Fermions. (c) Orbital valley 
splitting of LLs. (d) Additional Zeeman spin-splitting of LLs 
(only the n = level is shown.) 



Note, however, that the similarity exploited above with 
DF in QED is valid only in the vicinity of points Ki^2 
and is only approximate. In reality, the transformation 
properties of k with respect to the graphene crystallo- 
graphic group G (Table |T| are quite different from those 
of the real DF transforming with respect to the Lorentz 
group. This provides the additional contributions to 
QED-like terms, such as the triangular-wrapped nonlin- 
ear kinetic [27| and relativistic-noninvariant Coulomb in- 
teraction terms 
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The Bloch expansion ^ also remains valid in a mag- 
netic field H, although the slowly varying coefficients 
F{r) are now classified according to the discrete set of 
LLs (instead of continuous k). The usual way of intro- 
ducing jy = V X A consists of the Peierls substitution 



k 



k + -^A, 

c 



(9) 



which in the case of the Hamiltonian ^ is the same 
as replacing kx ± iky by the LL creation (annihilation) 
operators a^. 

The Hamiltonian has a set of discrete LLs having a 
square-root energy dependence on the level number n = 



3 



0, ±1, ±2, ... in a magnetic field (Fig[2]3): 



of the perturbation operator 



En = Seh\/'2v'^\e\hH;:\n\/c, Seh = sign(n), (10) 

which is quite different from the case of massive particles 
with En = hu!c{n + ^) > 0, {n = 0, 1, 2...), and has solu- 
tions with values both above and below the zero-energy 
LL, Eq = 0. The corresponding eigenfunctions can be 
written as an expansion over the n-th LL eigenfunctions 
fnir): 



F"-{r) = {cif\n\{r),Sehicif\n\-i{r), 

C2f\n\{r), SehiC2f\n\-l{r)}, 



(11) 



[note that f-i{r) = 0]. Each level, including n = 0, has 
the two-valley degeneracy provided by the complex con- 
stants ci,C2, the two-fold spin degeneracy and the infinite 
Landau degeneracy. 

Although the Peicrls substitution ([9]) conserves a rel- 
ativistic invariance of the Dirac equation in a magnetic 
field, the discrete crystal lattice background leads to an- 
other, weaker, requirement that the Hamiltonian of the 
system should be invariant with respect to the crystallo- 
graphic group of graphene in a magnetic field [28j : 



Gh — {CeR, C3, C2i?, cTxR, cTyR}. 



(12) 



In particular, this time non-invariant group Gh C G 
does not contain the translations Ti^2 that are incompat- 
ible with the translational magnetic group. The principal 
idea of the present work is that the graphene Hamilto- 
nian for charge carriers in magnetic field should have the 
more general form: 



H 



jI-IbH;, XfJ-BHz -va+ 
\ -"/hbHz -t)a" -XfiBHz ) 



(13) 



ifJ'B = \e\H/2mc ) that, besides the Peierls terms 
contains the "non-relativistic" A- and 7- corrections pro- 
vided by the orbital interaction of Bloch electrons with 
the magnetic field. 

These terms keep {^H"^) invariant under the opera- 
tion of group Gh and are produced by the matrix ele- 
ments: 



K2 



(14) 



IIJ'bHz 



and 



(15) 
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where m is the bare electron mass that accounts for the 
translational-invariant symmetry breakdown due to the 
discrete crystal background [29l . f30j . (Analogous terms 
for the time-symmetry-breaking field have been proposed 
in [31] for the orbital part of intrinsic spin-orbit coupling 
in graphene, which is minute.) 

The numerical parameter 7 of the matrix elements (jl4p 
is estimated in the tight-binding nearest neighbor ap- 
proximation (between sites A and B) as: 



7 ; 



t/2 



= 0.4, 



(17) 



where t — 3.033 eV is the tt-tt hopping integral, and 
a = 1A2A is the hexagon side (C-C interatomic dis- 
tance). It is more difficult to estimate the next nearest 
neighbor parameter A (between A and A') in (fT5|) . since 
the hopping integral falls off fairly slowly [t-^T^ (d) cx 1/ <P] 
but, clearly, A < 7. 

Diagonalization of the Hamiltonian can be easily 
done in terms of LL wave functions /„ (r) , presenting the 
resulting 4-component eigenfunctions F^{r) in the form 
(fTTj) with slightly different coefficients. This again gives 
the set of discrete LLs with n = 0,±1,±2,..., having 
the energies: 



En = Seh ^2vM^HMIc+{l±Xf ^ilHl (18) 
Special attention should be paid to zero LL n = with 
Eo = {-X±-i)iibH,, (19) 

and 

Fl{r) = {/o(r), 0, ±/o(r), 0}. (20) 

The new effect here, illustrated by Fig.[2j:, is the valley- 
splitting of each LL, marked by the ± sign and estimated 
as: 

Ai?„ ^ 7A f ff4) ' - 2 . lO-3^i/,[T]3/2 K, 

(21) 
(22) 



A£;o~27AiBiJ, ~1.37i/,[T]K. 



Being very small for non-zero LLs, n ^ 0, this splitting 
should be observable for zero LL in high fields. Note that 
this effect has purely orbital origin and is completely de- 
coupled from the additional Zeeman spin-splitting shown 
in Fig. Hi 

^E, ~ giiBH ~ 1.3|i/[T] K, 5 « 2, (23) 
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FIG. 3: Schematic of circular currents corresponding to Bloch 
functions *± = ± of ^Pli^ ^^ro Landau level. 



graphene. The effect occurs in a perturbative non-critical 
manner and is an intrinsic property of noninteracting 
fermions on a hexagonal lattice. As a consequence (ob- 
servable optically), the orbital splitting should not de- 
pend on the LL filling factor, unlike the result from other 
models. At the same time, the many-body and/or dis- 
order effects can amplify the orbital splitting (even for 
n 0), induce an additional symmetry breaking, and 
bring about a nontrivial field and filling factor depen- 
dence of the gap observed experimentally [ll|, ■ 

We grateful to Y. Kopelevich and M. Dyakonov for 
valuable discussions. IL thanks the support of the ANR 
agency (project LoMaCoQuP). 



because of very weak spin-orbital coupling. Unlike 
orbital-splitting, the spin-splitting is a function of the 
absolute value of H (rather than of Hz), that permits 
separating the two contributions AEg and AEn by their 
angular field dependence. 

To clarify the physical origin of the orbital splitting, 
consider the explicit form of the wave functions for car- 
riers located on the up/down shifted zero LL: 

$i(r) = (0^^,(0 = ^±{r)Mr), (24) 

i=l 

where the Bloch parts 

*±(0= *l(r)±vl>f^(r), (25) 

are presented in Fig. [3] and can be interpreted as a set 
of clockwise and counterclock-wise current loops circu- 
lating around every third hexagon. Therefore, it is the 
orbital paramagnetic interaction of these inter- valley cir- 
cular currents with that causes the splitting AEq. 

The current distribution for LLs with n ^ is more 
complicated since both clockwise and counterclock-wise 
current loops with different envelope LL functions [fn (r) 
and /„_i(r)] contribute to the wave functions of each 
split LL. The compensation of orbital momenta is almost 
complete, which explains the negligibly small splitting 
(I^Tj) of higher LLs. Note also that the additional contri- 
bution can be caused by the mtra- valley circular currents 
circulatin g ar ound next-nearest-neighbor plaquettes pro- 
posed in [l8|. Governed by the next-nearest neighbor 
parameter A, these currents do not contribute to n = 
LL splitting and contribute only very weakly to the split- 
ting of other LLs. Another consequence of the orbital LL 
splitting is the lattice period tripling produced by the 
network of circular current shown in Fig. [31 This field- 
induced breaking of graphene spatial symmetry can be 
observed for non-integer filling of zero LL, when clock- 
and counterclock-wise currents do not compensate each 
other. 

To conclude, we have proved that the orbital mech- 
anism is sufficient to explain the zero LL splitting in 
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